Multiple regions of quantum criticality in YbAgGe 

G. M. Schmiedeshoff, 1 E. D. Mun, 2 A. W. Lounsbury, 1 S. J. Tracy, 1 E. C. Palm, 3 
S. T. Hannahs, 3 J. -H. Park, 3 T. P. Murphy, 3 S. L. Bud'ko, 2 and P. C. Canfield 2 

1 Department of Physics, Occidental College, Los Angeles, California, 90041- 
mes Laboratory and Department of Physics and Astronomy, 
Iowa State University, Ames, Iowa, 50011. 
s National High Magnetic Field Laboratory, 
Florida State University, Tallahassee, Florida 32310. 
(Dated: January 13, 2013) 

Abstract 

Dilation and thermopower measurements on YbAgGe, a heavy-fermion antiferromagnet, clarify 
and refine the magnetic field-temperature (H-T) phase diagram and reveal a field-induced phase 
with T-linear resistivity. On the low-H side of this phase we find evidence for a first-order transition 
and suggest that YbAgGe at 4.5 T may be close to a quantum critical end point. On the high-if 
side our results are consistent with a second-order transition suppressed to a quantum critical point 
near 7.2 T. We discuss these results in light of global phase diagrams proposed for Kondo lattice 
systems. 

PACS numbers: 65.40.De,72.15.Jf,75.30.Kz,64.70.Tg 
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When a classical second-order phase transition is suppressed to zero temperature T by a 
tuning parameter (such as pressure, doping, or magnetic field H), a quantum critical point 
(QCP) can occur-; suppression of a first-order phase transition can lead to a quantum criti- 
cal end point (QCEP) 2 . In the vicinity of a QCP or QCEP quantum fluctuations associated 
with the zero-point energies of the adjacent T = phases can persist to remarkably high 
temperatures. These fluctuations can dramatically affect the interactions between parti- 
cles, leading to unusual thermodynamic and transport properties^"- and to novel states of 
matter- - -. 

YbAgGe, a stoichiometric heavy-fermion (HF) antiferromagnet, crystallizes in a hexag- 
onal ZrNiAl-type structure^. The zero-field electronic specific heat coefficient falls in the 
range 0.15-1.0 J/mol K 2 and the Kondo temperature is 20-25 K-. The Yb ions form a quasi- 
Kagome lattice in which magnetic coupling, geometric frustration and Kondo interactions 
compete in a manner which allows the suppression of low temperature AF order by modest 
applied magnetic fields, fields that tune the quantum critical behavior—^ and lead to a 
complex magnetic phase diagram^. 

The phase diagram of YbAgGe, summarizing earlier measurements^ with H perpen- 
dicular to the c-axis, is shown in Fig. [T^i. Solid (dashed) lines are guides to the eye for 
thermodynamic phase boundaries ('Hall lines' denoting features in the field- dependent Hall 
resistivity^ 1 ^) labeled by numerals 1-3 (4-6), lower case letters a-f label phases or regions 
of the phase diagram. Commensurate AF order is observed in the a-phase 1 - where a sharp 
first-order phase transition manifests along phase line 1-, and incommensurate AF order 
has been reported in the 6-phase^-. Relatively broad features manifest along the higher 
temperature part of phase line 3^, but below about 0.5 K the features sharpen and neutron 
scattering measurements reveal a return to commensurate AF order—. Based on the sup- 
pression of phase line 3 near 4.5 T and its near coincidence with Hall line 4 below 0.3 K, 
a field-induced QCP was proposed^^i. Intriguingly, several other features of this phase 
diagram remained mysterious. 

The low-T electrical resistivity is large^ and varies like T n with n ~ 1 in region d, smoothly 
increasing from ~ 1 to ~ 2 in region e, and ~ 2 in region f 19 . Such non-Fermi liquid (nFL) 
behavior is expected near a QCP- 1 ^, but the broad field range with n ~ 1 and the recovery 
of Fermi liquid (FL) behavior in fields so far above that of the QCP are surprising. Further, 
a logarithmic divergence of the specific heat appears most clearly for H ~ 7 T-, near Hall 
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FIG. 1: (color online) The phase diagram of YbAgGe below 1 K with H applied perpendicular to 
the c-axis. Panel (a) summarizes measurements from the literature, the variables carry their usual 
meanings. Panel (b) summarizes dilation and thermoelectric power data from this work with the 
same guides to the eye as panel (a) (see text). 

line 5, raising the possibility of (at least) one other phase being suppressed in fields well 
above that of the proposed QCP near 4.5 T. In this paper we describe longitudinal dilation 
and transverse thermoelectric power (TEP) measurements on YbAgGe that shed light on 
some of these mysteries and paint a fuller picture of quantum criticality in Yb-based HFs. 

Single crystals of YbAgGe were grown from an AgGe-rich ternary solution 9 . Longi- 
tudinal dilation and transverse TEP measurements were made with techniques described 
elsewhere^ 1 ^. H was applied normal to the c-axis. The coefficients of linear thermal ex- 
pansion and linear magnetostriction along the afr-axis (approximately parallel to [210]) are 
®ab = d(\nL ab )/dT and X a b = d{hiL ah ) / OH respectively where L ab is the thickness of the 
sample, 1-2 mm. The Seebeck coefficient is 5* = — AV/AT, where AV is the potential dif- 
ference across the sample when a temperature gradient is applied. Typical data are shown 
in Figs. |2] and [3] where nFL behavior is also illustrated. 
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FIG. 2: (color online) Representative longitudinal dilation data with H applied parallel to the 
aft-axis. Panel (a) shows the low temperature thermal expansion in the fixed fields listed next to 
arrows denoting the appropriate axes for each data set (the data at T/1.2 T have been shifted 
by +10/ — 9.0 x 10~ 6 /T for clarity). The Tj label features associated with the phase boundaries 
or crossovers of Fig. 1 (see text). A semi-log plot of a a b/T vs. T in the vicinity of 4.5 T is shown 
in the inset where the direction of increasing H is shown by a curved arrow and the solid line is a 
guide to the eye for the low temperature 4.7 T data. Panel (b) shows the magnetostriction at 387 
mK. The Hj label features in the data where j denotes the guides to the eye of Fig. la (see text). 
The inset shows an expanded view of the transitions below 6 T at 17 mK. 



In our dilation data, first-order phase transitions are identified by their peak-like shape 
and the presence of thermal or magnetic hysteresis (e.g. features labeled T a and i?i_4 in 
Fig. 2); second-order transitions do not show hysteresis within experimental uncertainty 
and usually show a step-like shape as illustrated by the dashed lines through the 3.4 and 
6.0 T data of Fig. 2a. Pronounced features in the 3.4 and 6.0 T data of Fig. 2a, labeled 
T C C 2 correlate well with extrema in X a b and are used to construct the phase diagram. The 
first-order transitions observed in the magnetostriction exhibit varying levels of magnetic 
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FIG. 3: (color online) Representative transverse TEP data with H applied parallel to the ab- 
axis and the temperature gradient in the basal plane, perpendicular to H. Panel (a) shows the 
temperature dependence of S across phase boundaries or crossover regions at the fields shown. The 
Tj label features associated with the phases (or regions) of Fig. 1. Panel (b) shows S/T vs. T in 
applied fields on a semi-log plot. The solid line through the 8 T data is a guide to the eye. The 
field dependence of S at 0.4 K is shown in the inset where the Hi label features in the data and i 
denotes the guides to the eye of Fig. la. 

hysteresis. This hysteresis becomes quite pronounced as the temperature falls below about 
0.2 K, as illustrated in the inset of Fig. 2b where the features associated with Hi and H% 
change dramatically in the dH/dt < data. For all of our T < 0.2 K magnetostriction data 
the sample was zero-field-cooled from temperatures above 0.3 K. Fig. lb was constructed 
from dH/dt > data. 

The phase diagram assembled from our dilation and TEP measurements is shown in Fig. 
lb where the solid and dashed lines are those of Fig. la, except for the dotted portion of 
phase line 3 where we observe no thermodynamic features. In comparison with the phase 
diagram of Fig. la we find: that phase line 2 joins phase line 3 to surround the c-phase; that 
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Hall lines 4-6 are all associated with thermodynamic phase transitions as T — > 0; and that 
phase transitions anchoring Hall lines 4 and 5 define the low temperature boundaries of a 
field-induced <i-phase (the high temperature limits of the n ~ 1 region of the resistivity agree 
reasonably well with the top of the (i-phase dome^) . The phase transitions are second-order 
across the tops of the c- and <i-phase domes and first-order along phase line 3 (Hall line 
4) below 0.3 K (0.2 K). The overlap of the phase line data determined from a a b, X a b, and 
S(T, H) is extensive along the c-phase boundary, there is some overlap on the high-field side 
of the (i-phase but no overlap is observed on the low-field side where a small gap appears 
between the a a b and X a b data, a gap Hall line 4 passes through. Broad extrema in a a b(T) 
and S(T), labeled with TJ, in Figs. and|3h. respectively, may represent a cross-over region 
between higher-T fluctuations and the lower-T incommensurate AF order of the 6-phase, 
leading to the uncertain nature of the phase diagram in this region. 

Several signatures of quantum criticality appear in our data. A change in the sign of the 
thermal expansion at a QCP has been predicted 2 ^, the phase diagram coordinates where 
a a b passes through zero are plotted as solid triangles in Fig. [lb. These 'lines of zeros' pass 
along the tops and high- if sides of the c- and <i-phases and presumably extend to T = 
near 4.5 and 7.2 T. Hall line 4 (5) correlates with the line of zeros associated with the c- 
phase ((i-phase) as T — > 0. A sign change in S(T) correlates with the field-induced QCP in 
YbRh 2 Si2 where it is attributed to an abrupt change in the Fermi surface 23 , though such a 
sign change does not appear to be a universal QCP signature^. S/T as T — » is predicted 
to reach its maximum value as the QCP is approached and its symmetry with respect to 
the QCP can help distinguish between theoretical models^. Sign changes in S(T, H) for 
YbAgGe are plotted as solid diamonds in Fig. lb and correlate with Hall line 4. Maxima 
in S(H) are plotted as solid stars in Fig. lb and correlate with Hall line 5 where the largest 
values of S/T as T — > also appear. 

The Griineisen parameter, characterizing the volume dependence of the energy scales in 
the system (which should be dominated by the quantum critical contribution as T — > 0), is 
defined as T a b = V m a a b/ kC p where V m is the molar volume 8 , k is the compressibility—, and 
C p is the specific heat 25 . T ab is predicted to diverge as a QCP is approached 26 , our results are 
shown in Fig. 4. The largest values of r a &, comparable to that of other HF compounds 27 , 
occur at 4.5 T. At this field the low-T upturn in T a b(T) is consistent with the onset of 
Griineisen divergence, perhaps with a power-law temperature dependence as suggested by 
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FIG. 4: (color online) The temperature dependence of the Griineisen parameter T a b of YbAgGe 
in applied magnetic fields. The curved arrow represents the direction of increasing field. (The 
apparent, slight downturn of the 8.0 T data at the low-T end is comparable to the noise in these 
data.) The inset shows the low temperature end of the 4.5 T data on a log-log plot. The solid line 
is a guide to the eye representing a T -5 / 2 temperature dependence. 

the data in the inset of Fig. 4 (though the temperature range over which divergent behavior 
is observed is too limited, at this time, for a definitive quantitative analysis). T a b(T) at 
7, 8, and 9 T also show low-T upturns suggesting that Griineisen divergence may develop 
at lower temperatures; we cannot, as yet, explain the large magnitudes nor the shapes of 
T a b(T) above 0.5 K in this field range. It will be necessary to extend these measurements to 
lower temperatures, higher fields, and at a higher density of field values. 

We take the presence of nFL behavior (both earlier— and current: the logarithmic tem- 
perature dependences of a a b/T and S/T suggested by the data in the inset of Fig. 2a and 
the main panel of Fig. 3b respectively), the features in the Hall resistivity^^ (Hall lines 4 
and 5), the zeros in a a b(T), and the features and zeros in S(T, H) as strong evidence for (at 
least) two regions of quantum criticality in YbAgGe near H& = 4.5 T (also supported by 
the onset of Griineisen divergence) and H c2 = 7.2 T. 

If a QCP at H& is due to the suppression of the AF transition characterized by the full 
phase line 3, as previously suggested^, then one would expect a continuous transition as 
T — > 1 , perhaps similar to that observed in YbRl^Sis 6 . However, the transitions along 
phase line 3 are first-order (see Fig. 2b) as T — > 0. A step-like feature developing in the low- 
T magnetization near H c i 25 and the thermodynamic structure of the c- and d-phases (two 
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second-order phase lines approaching/joining a first-order phase line extending to T = 0) 
are similar to that expected for the spin-flop class of metamagnetism^^. Metamagnets can 
exhibit a wide variety of multicritical behavior, including intermediate anisotropy scenarios 
leading to critical end points^. We suggest that YbAgGe may be close to a QCEP near if cl , 
possibly similar to that observed in S^R^O^ but with spin-flop metamagnetism driving 
the quantum criticality. 

The phase transitions at the base of Hall line 5 are clearly continuous as shown in the 
main panel of Fig. 2b while n ~ 1 and a a b(T) passes through zero here. In this region of the 
phase diagram both the specific heat and S/T are logarithmically divergent. We propose 
that a QCP occurs near H c2 = 7.2 T, driving the pronounced nFL behavior in this region 
of the phase diagram. 

Theoretical efforts, characterizing several quantum critical materials on a global phase 
diagram incorporating Kondo coupling and degree of magnetic frustration, suggest that 
YbAgGe may evolve from AF order (d-phase) through a spin-liquid phase (e-phase) before 
a FL (/-phase) is recoverecr^2 Recent work on other Yb-based HF compounds suggests^ 1 ^ 
that a quantum critical phase, bounded by QCPs at H C 2 and a hypothetical spin-liquid/heavy 
FL QCP— near 12 T, may underlie the spin-liquid. The linear T-dependence and large 
magnitude of the resistivity suggest strange-metal behavior for the c?-phase^, and as a spin- 
flop phase it will carry a net magnetic moment^, though the small step in M(H), about 
0.1/zb/Yd^, suggests that an underlying AF symmetry may still be present. The theoretical 
identification of the e-phase as a spin-liquid is supported by the large magnetostriction^ we 
observe. The proposed global phase diagram thus seems appropriate for YbAgGe. Elastic 
and inelastic neutron scattering would be profitable though challenging microscopic probes 
of these high- if phases. 

In conclusion, dilation and TEP measurements reveal high-if phase boundaries in 
YbAgGe that delineate the region of T-linear resistivity. On the low- if side this phase 
appears to be close to a QCEP near 4.5 T, associated with a first-order, most likely metam- 
agnetic, phase transition. On the high- if side this phase appears to end in the continuous 
suppression of a second-order transition ending in a QCP near 7.2 T, explaining the pro- 
nounced nFL behavior nearby. Even with the identification of this field-stabilized phase, 
there remains a clear nFL region over which the resistivity varies as T n with n continuously 
changing from ~ 1 to ~ 2 as H increases. Theory suggests a quantum critical phase and/or 
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spin-liquid in this region but more evidence, theoretical as well as experimental, is needed. 

This work was supported by the National Science Foundation under DMR-1006118. Work 
at Ames Laboratory was supported by the U.S. Department of Energy, Basic Energy Sci- 
ences, under Contract No. DE-AC02-07CH11358. Work at the National High Magnetic 
Field Laboratory was supported under the auspices of the National Science Foundation, the 
State of Florida, and the U.S. Department of Energy. 



1 S. Sachdev, Quantum Phase Transitions (Cambridge Univ. Press, New York, 1999). 

2 A. J. Millis, et al., Phys. Rev. Lett. 88, 217204 (2002). 

3 G. R. Stewart, Rev. Mod. Phys. 73, 797 (2001). 

4 H. v. Lohneysen, et al., Rev. Mod. Phys. 79, 1015 (2007). 

5 P. Coleman and A. J. Schofield, Nature 433, 226 (2005). 

6 P. Gegenwart, et al., Nat. Phys. 4, 186 (2008). 

7 A. J. Schofield, Phys. Stat. Sol. B 247, 563 (2010). 

8 R. P6ettgen, et al., Z. Kristallogr. 212, 58 (1997). 

9 S. L. Bud'ko, et al, Phys. Rev. B 69, 014415 (2004). 

10 Q. Si, Phys. St. Sol. B 247, 476 (2010). 

11 K.-S. Kim and C. Pepin, Phys. Rev. B 81, 205108 (2010). 

12 D. Gignoux and D. Schmitt, J. Alloys Compd. 326, 143 (2001). 

13 S. L. Bud'ko and P. C. Canfield, Physica B 403, 1230 (2008). 

14 S. L. Bud'ko, et al, Phys. Rev. B 71, 054408 (2005). 

15 S. L. Bud'ko, et al., Phys. Rev. B 72, 172413 (2005). 

16 B. Fak, et al., J. Phys. Cond. Matt. 17, 301 (2005). 

17 B. Fak, et al., Physica B 378-380, 669 (2006). 

18 D. F. McMorrow, et al., oral presentation, 25 th Int. Conf. on Low Temp. Phys., Amsterdam 
2008. 

19 P. G. Niklowitz, et al., Phys. Rev. B 73, 125101 (2006). 

20 G. M. Schmiedeshoff, et al., Rev. Sci. Instrum. 77, 123907 (2006). 

21 E. Mun, et al., Meas. Sci. Technol. 21, 055104 (2010). 

22 M. Garst and A. Rosch, Phys. Rev. B 72, 205129 (2005). 



9 



S. Hartmann, et al., Phys. Rev. Lett. 104, 096401 (2010). 

K. Sengupta, et al., Phys. Rev. B 81, 125129 (2010). 

Y. Tokiwa, et al., Phys. Rev. B 73, 094435 (2006). 

L. Zhu, et al., Phys. Rev. Lett. 91, 066404 (2003). 

A. de Visser, et al., Physica B 163, 49 (1990). 

E. Stryjewski and N. Giordano, Adv. Phys. 26, 487 (1977). 

A. Aharony, Lecture Notes in Physics 186, 209 (1983). 

I. Vilfan and S. Galam, Phys. Rev. B 34, 6428 (1986). 

S. A. Grigera, et al., Science 294, 329 (2001). 

J. Custers, et al, Phys. Rev. Lett. 104, 186402 (2010). 

P. Coleman and A. H. Nevidomskyy, J. Low Temp. Phys. 161, 182 (2010). 

J. McGreevy, Physics 3, 83 (2010). 

J. P. C. Ruff, et al, Phys. Rev. Lett. 105, 077203 (2010). 



10 



